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We perform a comparative study of applicability of the Multifractal Detrended Fluctuation Anal- 
ysis (MFDFA) and the Wavelet Transform Modulus Maxima (WTMM) method in proper detecting 
of mono- and multifractal character of data. We quantify the performance of both methods by using 
different sorts of artificial signals generated according to a few well-known exactly soluble mathe- 
matical models: monofractal fractional Brownian motion, bifractal Levy flights, and different sorts 
of multifractal binomial cascades. Our results show that in majority of situations in which one does 
not know a priori the fractal properties of a process, choosing MFDFA should be recommended. In 
particular, WTMM gives biased outcomes for the fractional Brownian motion with different values 
of Hurst exponent, indicating spurious multifractality. In some cases WTMM can also give different 
results if one applies different wavelets. We do not exclude using WTMM in real data analysis, 
but it occurs that while one may apply MFDFA in a more automatic fashion, WTMM has to be 
applied with care. In the second part of our work, we perform an analogous analysis on empirical 
data coming from the American and from the German stock market. For this data both methods 
detect rich multifractality in terms of broad f(ct), but MFDFA suggests that this multifractality is 
poorer than in the case of WTMM. 



I. INTRODUCTION 

It is well-known that the self-similarity of fractal struc- 
tures can be described by the so-called Holder exponents 
or the local Hurst exponents a. If the fractal is homoge- 
nous (monofractal) then it can be associated with only 
one Holder exponent, while in the case of a multifractal, 
different parts of the structure are characterized by dif- 
ferent values of a, leading to the existence of the whole 
spectrum f(a). In contrast to model fractals with a pre- 
cise scaling, many fractals, both the mathematical and 
the natural, reveal only a statistical scaling and this in 
particular refers to the so called fractal signals. An exem- 
plary temporal process with a trivial monofractal scaling 
is the fractional Gaussian noise; this process has only one 
a equal to 0.5 for the uncorrelated motion and a ^ 0.5 if 
any linear correlations exist in the signal. On the other 
hand, a process with either non-Gaussian fluctuations or 
with nonlinear temporal correlations can be multifrac- 
tal, and two or more values of a or even a continuous 
spectrum /(a) can be needed to characterize structure 
of such a process. 

In recent years much effort has been devoted to reliable 
identification of the multifractality in real data coming 
from such various fields like e.g. DNA sequences [10, OJ, 
physiology of human heart jj, |f| , neuron s piki ng |(| , at- 
mospheric science and climatology Ql la la 05 , financial 
markets JTl El El El El El 111 W El H, geo- 
physics [221 and many more. This is, however, a diffi- 
cult task mainly due to the fact that experimental data 
related to physiology, economy or climate is highly non- 
stationary and, additionally, the available data samples 
are usually rather small. One thus requires to apply 
methods which are insensitive to nonstationarities like 
trends and heteroskedasticity. In principle there are two 
competitive methods of detection of the multifractality 



which are commonly used in this context; both supposed 
to eliminate trends and concentrate on the analysis of 
fluctuations. Multifractal Detrended Fluctuation Anal- 
ysis (MFDFA) H IH |13 is based on the identification 
of scaling of the gth-order moments that power-law de- 
pend on the signal length and is a generalization of the 
standard DFA using only the second moment q = 2. 
The other method, Wavelet Transform Modulus Maxima 
(WTMM) [Hill consists in detection of scaling of 
the maxima lines of the continuous wavelet transform on 
different scales in the time-scale plane. This procedure 
is advocated as especially suitable for analyzing the non- 
stationary time series [27] . 

We therefore we shall test and compare the applicabil- 
ity of those two methods to the data coming from cither 
a few mathematical fractal models or empirical data col- 
lected from the stock market. In this work we represent 
a point of view of a practitioner who wants to choose a 
better tool for his analyses without entering subtle the- 
oretical considerations. A deeper analytical research, al- 
though important and desired in general, remains beyond 
the scope of this work. 

That data from various markets like the stock market, 
the foreign currency market and the commodity one are 
of multifractal nature , it is well-known from a numer- 
ous recent studies EllHE!ElE!ElElEiE!II|, 
that were carried out with the help of both methods. A 
few years ago the Multifractal Model of Asset Returns 
was developed [H |3(| in order to explain the origin 
of this multifractalit y. This model and its later modi- 
fications [371 l.'SSl l39l |40| use the multiplicative cascades 
which generate signals that are inherently multifractal 
and that are able to mimic some key properties of finan- 
cial data. The rationale behind the introduction of such 
a model was the observed correspondence between finan- 
cial market evolution and fluid turbulence 1411. From 
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the other point of view, the existence of the so-called 
financial stylized facts (fat tails of the fluctuation dis- 
tributions and lon g-la sting nonlinear correlations in the 
signals) (HiHilMlii 

can also be considered a source 
of the multifractal dynamics 0, 0, |2(J . 

Our paper is organized as follows: In Section 2, 
we briefly sketch the foundations of the MFDFA and 
WTMM methods. We apply them to a few types of 
model data like Brownian motion, Levy process and bi- 
nomial multiplicative cascades in Section 3. In Section 4 
we illustrate the performance of both methods in a con- 
text of real financial signals and, finally, we arrive at the 
concluding remarks in Section 5. 



The result of the MFDFA(7) procedure is the family of 
exponents h(q) (called the generalized Hurst exponents) 
which, for an actual multifractal signal, form a decreasing 
function of q, while for a monofractal h(q) = const. The 
singularity spectrum of the Holder exponents f(a) can 
easily be obtained from the generalized Hurst exponents 
by the following relations |47j 

a = h(q)+qh'(q) f(a)=q[a-h(q)) + l. (5) 

a characterizes the strength of singularities and /(a) can 
be considered the fractal dimension of a subset of the time 
series with singularities of strength equal to a. 



II. DESCRIPTION OF METHODS 



B. Wavelet Transform Modulus Maxima method 



A. Multifractal Detrended Fluctuation Analysis 

The Detrended Fluctuation Analysis [l| has recently 
become a commonly used tool in analyses of scal- 
ing properties of monofractal signals and in identify- 
ing correlations present in noisy nonstationary time se- 
ries [2^. The multifractal generalization of this proce- 
dure (MFDFA(O) [12 can be briefly sketched as follows. 
First, for a given time series cc(i), i = 1, N on a com- 
pact support, one calculates the integrated signal profile 
Y(j) 



Y(j) = J2(x(i)-<x>), j = l,. 



(1) 



i=i 



where < ... > denotes averaging over the time series, and 
then one divides it into M n segments of length n (n < N) 
starting from both the beginning and the end of the time 
series (i.e. 2M n such segments total). Each segment 
v has its own local trend that can be approximated by 
fitting an Zth order polynomial Pj and subtracted from 
the data; next, the variances for all the segments v and 
all segment lengths n have to be evaluated 



n 

F\u,n) = -Y,{Y[{v-l)n + J ]-Pi l \ ] )} 2 . 



(2) 



j'=i 



Finally, F 2 (v,n) is averaged over v's and the qih order 
fluctuation function is calculated for all possible segment 
lengths n: 



w = {^kE [F2 ^ n)]q/2 



2M„ 



t/<l 



gen. (3) 



The key property of F q (n) is that for a signal with fractal 
properties, it reveals power-law scaling within a signifi- 
cant range of n 



F q (n) - n 



h(q) 



(4) 



The alternative Wavelet Transform Modulus Maxima 
method is a technique based on the wavelet transform p4l 

MM 



)x(i) 



(0) 



where tp is a wavelet kernel shifted by n and s' is scale. 
The wavelet method can serve as a tool for decompos- 
ing the signal in time-scale plane; the resulting wavelet 
spectrum T^(n, s') can reveal a hierarchical structure of 
singularities (see Figure 1). As a criterion for the choice 
of the mother wavelet ip, a good localization in space 
and in frequency domains is recommended. The family 
of wavelets which is used most frequently in this case is 
the mth derivative of a Gaussian 



dx r 



(7) 



because it removes the signal trends that can be approx- 
imated by polynomials up to (m — l)th order |22| . 

In the presence of a singularity in data one observes 
the power law behaviour of the coefficients 



T^(n , S 'W Q ("°). 



(8) 



This relation, however, is not stable in the case of densely 
packed singularities; it is thus much better if one iden- 
tifies the local maxima of and then calculates the 
partition function from moduli of the maxima 

Z(q,s')= ]T |7V(n i «U')l 9 , (9) 
leL(s') 

where L(s') denotes the set of all maxima for scale s' and 
ni(s') stands for the position of a particular maximum. 
In order to preserve the monotonicity of Z(q, s') on s' , 
one has to impose an additional supremum condition 



Z(q,s')= ]T (supjT^nKs'OyOD* 
leL(s') s "- s ' 



(10) 
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For a signal with a fractal structure, we expect that 
Z(l> s ') ~ s' T ( q \ The singularity spectrum f(a) can 
now be obtained according to the following formulas 0] 

a = t' (q) and f(a) = qa — T(q). (11) 

Additionally, there is a relation between r(q) and the 
generalized Hurst exponents 

r(q) = qh(q) - 1. (12) 

Linear behaviour of r(q) indicates monofractality 
whereas nonlinear one suggests that a signal is multi- 
fractal. An example illustrating the WTMM procedure 
is shown in Figure 1 for the devil's staircase obtained by 
integrating the Cantor measure. It should be explained 
that we calculate a locally near each q and this can some- 
times lead to an unexpected shape of f(a) as it is seen 
in Figure 1(b). 

III. COMPUTER GENERATED DATA 

First we shall consider a few examples of signals asso- 
ciated with some well-known processes for which the ex- 
act theoretical results are available. Our objective is to 
confront the outcomes of the analyzed methods against 
theory and detect advantages or disadvantages of each 
of the two procedures. Our analysis of artificial data 
concentrates on several important issues regarding the 
performance of MFDFA and WTMM: 

• method's ability of correctly identifying monofrac- 
tal or multifractal character of signals 

• method's precision in evaluating /(a) spectra that 
agree with respective theoretical predictions 

• stability of results across different variants of 
MFDFA (different polynomials) and of WTMM 
(different wavelets) 

• stability of results across different realizations of a 
given stochastic process 

• method's ability of providing one with correct re- 
sults for short signals 

• quality of scaling in F q (s) and Z(q,s') and sensi- 
tivity of results to distinct choice of a fitting range 
of s or s'. 

For the sake of consistency, we did not apply WTMM 
to original signals but rather to their integrated versions; 
this allowed us to compare the results of WTMM with the 
ones of MFDFA which in fact also analyze the integrated 
signal profile Y(j) instead of x(j) (see Eq. Q). In order 
to be able to compare the results from both the F q and 
the Z function, we also have to derive 

\s'Z{q,s')\ 1 ' q ~s' h{ - q \ (13) 



It is noteworthy that for the data analyzed here we deal 
with the two main possible sources of multifractality: 
the nonlinear temporal correlations and the broad prob- 
ability density functions. Finally, we have to mention 
that an earlier approach to a comparison of MFDFA 
and WTMM |2j which was only briefly sketched there, 
showed that the former method can provide one with the 
results being in a better agreement with theoretical pre- 
dictions than does the latter one. However, such a com- 
parison was not a central issue of the cited work which 
was rather concentrated on the performance of MFDFA 
itself. Here we present a more thorough study on this 
subject. 



A. Monofractal signals: Brownian motion 

We start our comparative study with testing each 
method's ability of identifying and quantifying monofrac- 
tal data. This is of great importance since in many prac- 
tical situations the question which is usually addressed 
first is whether data under study is mono- or multifrac- 
tal. In general, even if a method works satisfactorily well 
for very long signals, it may happen that for typical ex- 
perimental time series the finite-size effects would lead to 
broadening of the f(a) spectra which are no longer point- 
like and this, obviously, may cause a spurious detection of 
multifractality. We would prefer a method which works 
more robust in such situations. 



1. Classical Brownian motion 

The first data type which we investigate is a simple case 
of a monofractal time series represented by the classical 
Brownian motion with the Hurst exponent H = 0.5. This 
process can be classified as stochastic with the stationary, 
independent and Gaussian-distributed increments. In or- 
der to obtain statistically significant results we carry out 
calculations on K = 10 independent realizations of this 
process. Theoretical spectrum consists here of one point 
localized at a = 0.5 and /(a) = 1 (Eq.©). 

Uppermost panel of Figure 2 shows exemplary plots of 
the fluctuation function (Eq. QJ, open symbols) and the 
rescaled partition function (Eq. (|13fl . filled symbols) for 
a single realization of the process and for a few different 
values of the Renyi parameter q. As it can be seen, high 
absolute values of q typically correspond to plots with 
worse scaling, while small |g|'s are associated with plots 
exhibiting more clear power-law shape. This observation 
is valid for both MFDFA and WTMM, but in this par- 
ticular case of a Brownian process, scaling for q <C is 
worser in WTMM than it is in MFDFA. This may cause 
some ambiguity in choosing a plot range for which the 
fitting procedure according to Eq. I|13fl is applied. This 
ambiguity can affect determining r(q) and thus can in- 
troduce undesired values of a ^ 0.5. 
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Taking all independent process realizations into our 
consideration, we derive the mean multifractal spectra 

K 

^ = Y j r {k \q), (14) 
fc=l 

where the average is taken over all individual data sam- 
ples. The so-calculated spectra (symbols) can be seen 
in bottom part of Figure 2 together with the theoretical 
linear spectrum (solid line) . We applied the third deriva- 
tive of a Gaussian (?/> 3 ) which is orthogonal, thus insen- 
sitive, to quadratic trends in a signal. For consistency, 
in MFDFA we chose polynomials P (2) (MFDFA(2)) and 
therefore we were also able to remove trends up to 
quadratic one. 

An issue which has to be examined in this context is 
the range of q used in the analysis. On the one hand, 
it should be as wide as possible in order for the method 
to be capable of detecting even subtle multifractal ef- 
fects in a signal. Thus, by taking an extremely narrow 
range of the Renyi parameter one might see a false, al- 
most point-like f(a) even for a signal actually compris- 
ing a variety of singularities of different strength. On the 
other hand, taking a too large \q\ can produce statisti- 
cally meaningless results based on an insufficient number 
of time series points. Figure 3 shows how the mean sin- 
gularity spectra f(a) calculated from r(q) according to 
formula (|11|) for the classical Brownian motion depend 
on (jr-range used: — 3 < q < 3 (top), — 5 < q < 5 (mid- 
dle) and — 10 < q < 10 (bottom). The most striking 
feature is a broad f(a) calculated by means of WTMM 
(right column) as compared to its MFDFA counterpart 
(left column) . It is interesting that only MFDFA gives a 
perfect identification of monofractality for the most nar- 
row interval of \q\ < 3, while WTMM does not. In the 
middle panels we see that while MFDFA still offers the 
spectrum that can approximately be considered a sin- 
gle point, WTMM shows f(a) of a significant non-zero 
width Aa := a(q m - nl ) — a(q max ) that can even suggest a 
sort of multifractal scaling. By extending q up to ±10, 
WTMM completely fails: the shape of the corresponding 
singularity spectrum becomes evidently multifractal. In 
contrast, despite its parabollic shape, f(a) for MFDFA 
is still defined on a relatively narrow Aa not allowing 
us to erroneously assume a multifractal character of the 
u nderly ing process. In each panel, standard deviations 
of /(a) are denoted by horizontal and vertical error bars. 

There is no simple and straightforward explanation of 
how such a discrepancy between MFDFA and WTMM 
does occur. We have already seen in Figure 2 some prob- 
lems with scaling of the partition function Z(q, s') for 
g's distant from zero. However, this lack of an ideal 
power-law behaviour cannot fully account for the ob- 
served strong deviation of /(a) from its theoretically 
predicted form. This is because we carried out analo- 
gous calculations for a few different possible intervals of 
scales in which the fitting of Z(q, s') was done according 
to Ea. l|13[l . and in each case we eventually arrived at the 



qualitatively similar /(a) spectra. Another source of er- 
ror can potentially be a poor statistics of data. Let us 
rewrite Eq. (|10|) in the following form 

Z(q,s')= £ 6( S ') 9 (15) 
leL(s') 

and then let us define an effective number of maxima 
used in the calculation of Z(q, s') for a fixed q and s' and 
express it by a fraction of the total number of detected 
maxima: 

i 

> 0.9- Z(q,s% (16) 

where / and L(s') have the same meaning as in Eq. l|10|) 
and A = #L(s'). We also define a similar quantity for 
MFDFA: 

R[jl a (n) = — min Mv : (— V \F 2 (n, v)] q / 2 ) 1 ^ 

> 0.9 -F q (n)}, (17) 

where v, M n and F 2 are the same as in Eq.©. Both 
these quantities are presented in Figure 4 for three 
choices of q < and for a few different process real- 
izations. As the plots show, for q = —3 and q = —5 
both the estimates of Z(q, s') and of F q (n) are based on 
a significant fraction of data for a vast majority of s' or 
n values. In contrast, for q = — 10 the statistics is evi- 
dently poorer and might be considered unsatisfactory. In 
this case a too small fraction of the number of maxima 
(even < 10 maxima in absolute numbers) contributes to 
the scaling functions F q (n) and Z(q, s') and thus they de- 
scribe only the scaling properties of a few largest events 
instead of the properties of the whole signals. Curiously, 
despite the fact that both A and RS$)„, ar e sta- 

" MFDFA WTMM 

tistically significant for q — — 5 and q = — 3, the output 
of MFDFA and WTMM is completely different. Consis- 
tently, qualitatively similar findings were also obtained 
for the other types of monofractal processes studied in 
this work. Thus we are justified to conclude that data 
statistics cannot be among principal sources of discrep- 
ancy in precision of both methods. 

Based on these outcomes, from now on we will restrict 
our analysis to — 5 < q < 5. We sample this parameter 
with Aq = 0.1 frequency for \q\ < 3 and with Aq = 0.5 
for \q\ > 3. 

The above discussion might indicate that a problem 
with the incorrect quantifying monofractal processes by 
WTMM as compared to MFDFA can rather be inherent 
to the WTMM procedure. For example, it may require 
significantly longer time series than does the detrended 
fluctuation method in order to obtain a good convergence 
in Eq. (|10l) for larger |g|'s. Therefore, in order to investi- 
gate how the methods work in respect to sample size, 
for the same type of Brownian motion we create sets 
of time series of different lengths (N — 15, 000, N — 
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65,000, N = 130,000 data points). As Figure 5 shows, 
we observe a noticeable ^-dependence of /(a) especially 
for WTMM: the longer the time series, the better agree- 
ment with theory. It can also be seen that for smaller N 
the spectra produced by WTMM are much more unsta- 
ble in terms of standard deviations than are the spectra 
for MFDFA. Obviously, standard deviations reveal also 
a strong dependence on N. Thus, for e.g. N = 15, 000, 
/(a) for a single process realization is likely to falsely 
indicate multifractality and temporal correlations (per- 
sistence or antipersistence) in a completely random, un- 
correlated signal. MFDFA seems to be more powerful 
here since even for short signals the standard deviations 
are acceptably small and the maximum of the average 
spectrum lies almost ideally at a — 0.5. In contrast, for 
the signals as long as N = 130, 000, WTMM provides us 
with the results that are not satisfactory. Perhaps, for 
much longer signals of length exceeding 10 6 data points, 
the situation would improve but, from a practical point 
of view, empirical signals of such a length are rarely avail- 
able. This strongly favours MFDFA as a monofractality 
detector. 

Next, we expect the results of an analysis to be insen- 
sitive as much as possible to a choice of the detrending 
polynomial in MFDFA and of the Gaussian deriva- 
tives i/A™* 1 in WTMM. Any significant dependence would 
limit robustness of a method since in the case of a real 
financial market one does not know a priori which poly- 
nomial or which wavelet function can be optimal. Fig- 
ure 6 shows the mean singularity spectra /(a) obtained 
by using MFDFA(Z) with I = 1,2,3,4. For pW we do 
not observe exactly single-point spectrum but rather a 
very narrow parabola. By increasing / we see decreasing 
Aa, the effect that is easily understandable: higher-order 
polynomials can better detrend the data. It seems that 
even P^ works sufficiently well in the present case. Suc- 
cessfully, in all cases the location of the spectra agrees 
perfectly with the expected H = 0.5. The correspond- 
ing results for WTMM (m = 1,2,3,4) are collected in 
Figure 7. We see that although maxima of f(a) are sit- 
uated correctly, the spectra are definitely too wide to 
be considered monofractal. The largest discrepancy be- 
tween these results and the theory is for tfj^, while for 
the higher-order derivatives the spectra are comparable 
in their widths (although still far from being monofrac- 
tal). It is also noteworthy that starting from m = 2, 
the standard deviations tend to increase with increasing 
derivative order. For MFDFA, where the standard devia- 
tions are much smaller, such a behaviour is not observed. 



2. Fractional Brownian motion 

In contrast to the above uncorrelated ordinary Brown- 
ian motion, the fractional Brownian motion with H ^ 0.5 
can serve as an example of a monofractal process with 
temporal correlations. Impact of these correlations on 
the results of MFDFA and WTMM can be evaluated from 



Figures 8-10. The two following cases are discussed: an 
antipersistent process with H = 0.3 and a persistent one 
with H = 0.75. Comparison of the /(a) spectra in Fig- 
ures 8 and 9 calculated for different PW (MFDFA, left 
columns) and ij)(m) (WTMM, right columns) leads to 
esentially similar conclusions as in the case of H = 0.5. 

(a) In each case MFDFA acts excellently for the an- 
tipersistent process: it shows clearly a complete lack of 
multiscaling. In the persistent case the spectra are less 
perfect, although they still cannot be erroneously con- 
sidered multifractal. As compared to MFDFA, Aa for 
WTMM is large and thus the corresponding WTMM 
spectra fail to reflect the actual monofractality of the 
signals. However, for high Gaussian derivatives like 
WTMM spectrum can be very narrow and its Aa can 
potentially be used as an indicator of possible monofrac- 
tality of an analyzed signal. 

(b) MFDFA provides us with the correct value of H 
in both the persistent and the antipersistent case. From 
this point of view WTMM also works well as regards the 
average spectra, but it turns less reliable for the individ- 
ual signals. 

(c) The average outcomes of MFDFA can be accounted 
stable across their different-order variants, while WTMM 
tends to work better if one uses higher Gaussian deriva- 
tives (to > 3). 

It should be stressed that we assume no a priori knowl- 
edge of the processes underlying the data. This, for ex- 
ample, distinguishes our analysis from the one carried out 
by authors of ref. who also applied WTMM to the 
fractional Brownian motion with H = 0.3. They took 
an advantage of the knowledge that the analyzed pro- 
cess had to reveal monofractal scaling and, accordingly, 
they were able to fit a straight line to the computed r(q). 
As a consequence, they obtained the exact single-point 
/(a). However, for data with unknown properties one 
has to calculate spectra by local fits and for WTMM this 
leads to a dispersion of a's even for monofractal pro- 
cesses. On the other hand, authors of ref. [22] carried 
out an MFDFA analysis based on three types of Brow- 
nian motion (H = 0.25, H = 0.5, H = 0.75) and they 
obtained similar results to the ones presented here. Even 
more, they were able to show that MFDFA can be reli- 
able even for \q\ > 10. 

Finally, we take a look at the results obtained for dif- 
ferent time series lengths (Figure 10). For the compari- 
son purpose we choose I = 2 and m = 3. For H = 0.3 
we see that while MFDFA gives us a narrow spectrum 
even for a signal as short as N = 15, 000 with increas- 
ing its accuracy with increasing N, the WTMM method 
occurs totally unreliable for shorter signals. The case of 
H = 0.75 does not differ much from the previous one but 
now MFDFA slightly looses its perfect accuracy and does 
not converge so well to a monofractal spectrum with in- 
creasing N. This is especially evident if a broader range 
| q | < 10 is used (not shown). Broadening of the MFDFA 
spectrum is a gradual process that becomes more evident 
for persistent signals with H 3> 0.5. It is worth mention- 
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ing t hat for both methods and for all N the maxima of 
f(a) point at the correct values of a. 

B. Bifractal signals: Levy processes 

The second interesting class of processes which are of 
high practical utility are the stable Levy processes. Their 
applications range from physiology to financial markets; 
in the latter case it is hypothesized that the modification 
of this kind of processes, i.e. the truncated Levy flights 
(with the exponent-suppressed Levy distribution tails) 
describes the price fluctuations of stocks. The scaling 
exponent r(q) for the truncated Levy flights (with the 
Levy parameter ol) and for non-Levy signals with the 
tails obeying the power law distribution P(x) ~ x -( a L+i) 
can be expressed by |U H3 

r(q) = l q/aL ~ 1 ( ( q - aL \ (18) 
W \ (q>a L ) y ' 

and the associated singularity spectrum by 

ol = ( X/ " L { ( q J UL \ f(a) = ( I <? J aL \ . (19) 
\ (q > a L ) Jy ' \ (q> a L ) y ' 

As these expressions show, signals with the truncated 
Levy distributions are rather bifractal than multifrac- 
tal |28j . From a dynamical point of view, the bifractal 
spectra or, more generally, the spectra revealing singu- 
larity in r(q), can often be seen in systems exhibiting 
phase transitions like e.g. chaotic systems with intermit- 
tency |2^, |3(], El H3 ■ It should also be noted that in 
principle the moments higher than do not exist at all, 
but as we consider the time series of finite lengths, we can 
also calculate h(q) for q > aL- We chose a heavy-tailed 
distribution with ar, — 1.5 and generated K = 10 time 
series of length N = 250, 000 data point s eac h. Examples 
of F q (s) and \s'Z(q, s')| 1/<? as well as of r(q) for MFDFA 
(/ = 2) and for WTMM (m = 3) are shown in Figure 
11. Scaling of F q (s) worsens for large positive q's (tri- 
angles down), while scaling of \s'Z(q, s')] 1 ^ proves weak 
for strongly negative q J s. 

Figure 12 presents results for different variants of 
MFDFA and WTMM. It is interesting to note that 
MFDFA (left column) presents a good agreement with 
theory near the bifractal points fEa. H19|) ) where there 
are two dense clusters of symbols. Unexpectedly, for 
each P(0 we also observe a continuous transition between 
these clusters. However, such a spurious transition is 
rather inevitable for real data supposed to exhibit two 
different linear regimes of r(q) (see also j2j|[3(j for other 
examples of spectra displaying an analogous transition). 
Main difference between spectra for different values of 
I is that the cluster at a = 1/or, tends to be better 
localized for I > 1. On the other hand, for all Gaus- 
sian derivatives in WTMM we observe roughly the same 
artificial transition as for MFDFA, but there is also an 
additional spurious arm of the spectrum for large a (i.e. 



strongly negative q) that effectively disperses the spec- 
trum beyond its theoretical prediction. This arm can be 
a consequence of weak scaling in Figure 11 that in turn 
can be a consequence of inherent problems with correct 
calculation of the multifractal spectra for q < 0. In this 
case estimation of F q (n) or Z(q, s') might be strongly bi- 
ased by the existence of very small values in time series 
which can be largely amplified if their negative exponents 
are considered; thus they can dominate the results com- 
pletely Unfortunately, this situation only slowly 
improves with a time series length so that some other, 
more refined techniques are required [34|. 

Influence of N on the results can be inferred after in- 
specting Figure 13. Apart from the longer error bars for 
N = 100, 000, for both methods the spectra for longer 
signals are smoother near the theoretical points, which 
reflects more reliable fits according to Eq. @ and Eq. (|13fl . 

C. Multifractal signals: binomial cascades 

In this subsection we present results for processes 
which are inherently multifractal, i.e. binomial mul- 
tiplicative cascades [22, 0- Processes of this kind 
are commonly used to model fluid turbulence and due 
to the recently-formulated hypothesis of similarity be- 
tween the turbulence and the evolution of financial mar- 
kets |4l) , th ey are more and more often applied in econo- 
physics [35l l36l l38l l39l l40| . 

1. Deterministic case 

Let us consider a probability measure /xq distributed 
uniformly on interval [0, 1] and two numbers mo, mi such 
that mo + mi = 1. In the first step we uniformly spread 
a fraction mo of total mass on the left subinterval [0, 1/2] 
and a fraction mi on the right subinterval [1/2, 1]. In the 
second and subsequent steps we repeat this procedure for 
each of the subintervals using in each step the same frac- 
tions mo and mi of the higher-level subinterval mass. In 
fcth step, each subinterval j (j = 1, 2 k ) can be la- 
belled by a unique sequence = T)i^r$f' ' ...r]^\ where 
■q\ ' is either or 1. Thus, a measure of this subinter- 
val is /ife^^'] = mj ^m™ , where n(j) denotes 
a number of unities in the binary representation of j. 
This construction leads to preservation of the total mass 
on subinterval [0,1], i.e. the measure fik is conservative 
fJ>k[v ] = 1- I n the limit k — > 00 the measure /i^ 
goes to binomial measure /x. 

The above-defined procedure generates a binomial cas- 
cade that after fc max steps can be represented by a time 
series {xj}jLi of length N — 2 femax such that 

Xj = a n( > j - x \l - a) fe — -««-!), (20) 

where a = mo and a G (0.5, 1). The resulting signal pos- 
sesses singularities of strength depending on the param- 
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eter a and, for a significantly less than 1, its multifrac- 
tality comes mostly from the temporal correlations (for 
a — ► 1 the broad probability distribution of Xj also con- 
tributes much) . The analytical expression for the scaling 
exponent and for the singularity spectrum can both be 
derived straightforwardly |22| : 



r(q) 



ln(2) 



•ln[o 9 + (1 -af 
M2) 



ln(a) + (l-a) g ln(l-q) 
a<? + (1 - a)i 



(21) 



(22) 



/(«) = - 



q a? ln(a) + (1 - a) q ln(l - a) 
ln(2) a? + (1 - a)? 

-ln[a 9 + (1 - a) 9 ] 

M2) ■ 



(23) 



In order to check how the MFDFA and WTMM meth- 
ods work for different values of the parameter a, we con- 
sider time series constructed for a = 0.55 (a relatively 
smooth signal), for a = 0.75 (existence of sharp singular- 
ities) and for an intermediate case of a = 0.65. It comes 
from the definition (Eq. i|2L)fll that by increasing a we en- 
hance the role of heavy-tailed p.d.f. of {Xj} and we also 
make the multifractality richer (larger Aa). Conversely, 
for a — > 0.5 the theoretical singularity spectrum tends 
to a monofractal one. We perform calculations on time 
series of length N = 131,072, i.e. in each case we stop 
the cascade-generating procedure at fc max = 17. Unlike 
for the Brownian and Levy processes, here due to the de- 
terministic nature of the cascades under study we create 
only one time series for each value of the parameter a. 
Figure 14(a) shows exemplary plots of F g (s) for MFDFA 
(PW) and |s'Z( g ,s')| 1/9 for WTMM (^ (3) ). Globally for 
a large range of s the scaling is relatively good despite 
the local fluctuations. Similar results can be obtained by 
using different polynomials in MFDFA or different Gaus- 
sian derivatives in WTMM (not shown here) . As regards 
r(q) shown in the panels (b) and (c), the estimated spec- 
tra resemble the theoretical one with the WTMM-based 
spectrum presenting almost perfect agreement with the- 
ory. 

The singularity spectra computed for distinct variants 
of MFDFA and WTMM are compared in Figure 15 with 
the ones theoretically derived according to Eq. (|22H and 
Eq. (|23|l . Interestingly, for a = 0.55 the only case in 
which we observe an agreement between theory and nu- 
merical estimates is for WTMM with ip^> (top right 
panel). For other wavelets the spectra deviate in either 
a or in /(a) and the same happens for all the analyzed 
polynomials of MFDFA (top left). A much better per- 
formance of WTMM is seen for a = 0.65 (middle right) 
where only %p^ does not reproduce the theoretical spec- 
trum, while by using the other Gaussian derivatives we 
get the correct outcomes. MFDFA offers us the f(a) 
spectra that slightly deviate from the expected one for 



positive values of the Renyi parameter q (corresponding 
to large signal fluctuations, left arms in middle-left panel 
of Figure 15) but their shape does not depend on the 
polynomial . Similar conclusions regarding MFDFA 
can be formulated for a — 0.75 (bottom left). In contrast, 
WTMM seems to be unable to cope with the latter kind 
of data for small (negative) q's (small signal fluctuations) : 
it significantly underestimates width of the spectrum on 
higher-a side which is most strikingly evident for tp^ 
and ip^ (bottom right). An agreement between theory 
and practice is reached only for positive q's except the 
wavelet ip^ . In general, as regards the whole spectra, 
for a — 0.75 the wavelets tp^ 1 ' an d V 1 behave more 
reliably than the remaining ones. It is also noteworthy 
that for all the studied values of a, the MFDFA results 
are remarkably independent of the polynomial choice, the 
property that cannot be attributed to WTMM and the 
choice of a wavelet. 

For next point of our analysis we prepare time se- 
ries of different lengths N = 16, 384, N = 65, 536 and 
N = 131, 072 in order to investigate how smaller N can 
affect the reliability of calculations as compared to long 
signals. Based on Figure 15, for each a we choose a 
wavelet with the best performance: tp^ f° r a = 0.55, 
^(3) for a = 0.65 and V (2) for a = 0.75. MFDFA is 
represented by the polynomial p( 2 \ Figure 16 collects 
plots of the corresponding singularity spectra; in a sharp 
contrast to Figures 10 and 13, in the present case the 
spectra are practically not sensitive to time series length 
in a wide range of N. At the end, we note that the out- 
comes presented here for MFDFA agree with the ones 
from reference 221. 



2. Stochastic case 

The above-described procedure can be generalized in a 
number of ways, e.g. by increasing the number of equal- 
size subintervals into which an interval is splitted at each 
cascade stage k and/or by randomizing the mass alloca- 
tion among the subintervals. This latter way is especially 
appealing because it allows one to eliminate the unreal- 
istic determinism of the previous example. Even more 
general models that can be related to financial data are 
available in literature (see e.g. H3,ESE3l) but within the 
scope of the present work it is sufficient that we concen- 
trate exclusively on some simplest stochastic cascade pro- 
cesses of multifractal nature. Proceeding along this line, 
we preserve the binomial character of the cascade but we 
replace the conservative measure with a multiplicative 
one that allows the multipliers rrii to be independent, 
identically distributed random variables drawn from a 
specific distribution. Now the total mass on interval [0, 1] 
is preserved only in a statistical sense: E(J2i mj) = 1. It 
has been shown that a resulting OTOcessproduces sig- 
nal with multifractal properties |36l l50l l5lj . Random 
multiplicative cascades, and especially their iterative ver- 
sions [39t |52j , are in practice much more interesting than 
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their deterministic counterparts due to the fact that they 
are able to mimic stochastic character of financial volatil- 
ity fluctuations over time. Here we discuss the singularity 
spectra of three common versions of such processes, i.e. 
Poisson, Gaussian and gamma multiplicative cascades. 

a. Log-Poisson cascade Let us start with a discrete 
cascade characterized by the multipliers m(r)i, r]k) 
such that M(r]i, ...,r)k) ■— — log m(rji, ...,r)k) have Pois- 
son distribution 



p(x) 



e 7 7 a: 



At a stage k, an interval j has the mass 

Mfo^D = m (Vi )m (vi, m) ■ ■■■ ■ m (Vi, --,Vk) 



and therefore 



ln// fe = ^M(r?i 



(24) 



(25) 



(26) 



The sum on the right-hand side is also Poisson- 
distributed with 7 — > kj in Eq. (|24|) and this leads to 
the following formula for the singularity spectrum |36l | 



/(«) 



1 



7 
In 2 



alog 2 (je/a). 



(27) 



This spectrum assumes its maximum at ao = 7 and 
has negative values for a ^ 7 and for a if 
7 > In 2. Based on the above-described generating pro- 
cedure we create a time series representing the k = 17th 
stage of the log-Poisson binomial cascade and by apply- 
ing MFDFA and WTMM we estimate f(a). Figure 17 
collects the results obtained with different variants of 
wavelets (ip^K ?/^ 4 ' ) ) and polynomials {P^\ P^). 

For WTMM (right-hand side) we observe the mean 
spectra whose increasing arms generally agree with the 
spectrum of Eq. 1)2 7fl but, on the contrary, the decreas- 
ing arms completely fail to comply with theory. MFDFA 
works substantially better than WTMM for a > ao, even 
if it does not reproduce the theoretical curve ideally after 
averaging f(a) over all 10 individual cascade realizations. 
This resembles the results for the deterministic cascades 
with a = 0.75 (Figure 15), where we also noticed prob- 
lems with estimating f(a) for a > 1 by means of WTMM. 
Therefore, this seems to be a more general issue of deal- 
ing with signals comprising relatively smooth singulari- 
ties associated with high a's. Both for MFDFA and for 
WTMM, significant error bars on both coordinates reveal 
strong instability of the calculated spectr a whi ch vary 
from sample to sample. Figure 18 displays /(a) for dif- 
ferent time series lengths (different k's). Again, MFDFA 
gives better results in respect to Eq. lt77|) even for rela- 
tively short signals than does the wavelet-based method. 
These outcomes for the log-Poisson cascades collected in 
Figures 17 and 18 qualitatively resemble the results of 
ref. for the same type of data. 



b. Log-gamma cascade As another example of a 
multiplicative cascade we investigate a process with the 
multipliers whose logarithms are taken (with minus sign) 
from the gamma p.d.f. 



p(x) = (Fx 



/r( 7 ) 



(28) 



where /3, 7 > 0. A sum of k i.i.d. random variables 
with such distributions is also a gamma distribution with 
7 — » try. Following |36| we can write the analytical form 
of the /(a) spectrum 



f(a) = l + 7 log 2 (a/3/7) + (7-a/3)/ln2 



(29) 



which reaches its maximum for ao =7//?. 

The plots in Figures 19 and 20 qualitatively resemble 
those from Figure 17 and 18 and, at least in part, similar 
conclusions can be drawn in the present case as it was 
for the log-Poisson one. Though, we note that at the 
decreasing part of the /(a) curves the error bars indicate 
extreme unstability of the computed spectra, suggesting 
that here neither of the two methods can be trusted for 
q < 0. This is caused by the fact that in this case plots of 
F q (n) and of \s'Z(q, s')\ 1 / q (not shown) present so weak 
scaling that even a small change of the fitting range of n 
and s' profoundly affects f(a). (See also the respective 
discussion related to the Levy processes.) 

c. Log-normal cascade Finally, we consider the mul- 
tipliers whose logarithms are taken from a Gaussian dis- 
tribution 



p{x) = 



1 



(2ira 2 



The singularity spectrum for a log-normal cascade 



/(«) = 1 - 



1 



2 In 2a 2 



(a - A) 2 



(30) 



(31) 



has a maximum for ao — X and extends over — 00 < a < 
+00 due to a possible lack of mass conservation. This 
theoretical spectrum has to be compared with spectra 
derived from time series by means of the MFDFA and 
WTMM methods. Figure 21 present such a comparison 
for the polynomials f , ' 1 ),..,,pW used in different vari- 
ants of MFDFA (left column) and for a few Gaussian 
derivatives ■..,ip < " 4 ' > used in WTMM (right column). 
It is clear that, as regards the mean spectra calculated 
from 10 independent realizations of the cascade and a 
range of a's for which /(a) > 0, both methods give re- 
sults that for a log-normal cascade are much closer to 
theoretical predictions than it was for the two cascades 
discussed above. This is particularly evident for the de- 
creasing part of the spectra, despite the existing devia- 
tions for a > 1.5. We observe a satisfying stability of the 
mean spectra across different variants of the methods and 
we do not notice any significant differences between the 
spectra produced by MFDFA and by WTMM except the 
fact that the standard deviations of data points tend to 
be larger for WTMM than for MFDFA. Also if we look at 
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the average spectra evaluated for different signal lengths 
(Figure 22) we see that the results differ mainly in stabil- 
ity across process realizations. In general, the standard 
deviations in a are relatively small while they are large in 
/(a); they are also much smaller in the increasing parts 
of the spectra than in the decreasing ones. 



IV. STOCK MARKET DATA 

In this Section we apply both methods of the multi- 
fractal analysis to real data from a stock market. Let us 
denote by P S (U) the price of an asset s at the ith con- 
secutive time instant (i, may in general not be equally 
spaced in time). We create a time series of logarith- 
mic price increments = ln(P s (i i+1 )) — ln(P s (£,)), 
where i = I, ...,N. The difference between our price in- 
crements p s (U) and the returns is that ti denote the mo- 
ments of transactions instead of the moments of constant- 
frequency data sampling. It should be noted that such 
a choice of data, motivated by our recent interest [20|. 
leads to the qualitative results that can be also drawn 
for the standard returns. 

We carried out the calculations for the time series 
representing 30 stocks comprised by Deutsche Aktienin- 
dex (DAX) and for the ones representing 30 Dow Jones 
Industrial stocks (DJIs) in the period 12/01/1997 - 
12/31/1999 53]. Each time scries under study had differ- 
ent length ranging from 63,000 (Karstadt) to 2,236,000 
(The Walt Disney Co.) data points. In a preliminary 
stage the data was preprocessed in the following way: we 
removed all overnight price increments (they were char- 
acterized by different statistical properties than the rest 
of the increments) and we also removed all the constant 
price intervals with more than 20 consecutive zeros (the 
MFDFA procedure requires a compact support of the 
signal in order to avoid the divergence of F q for q < 
and, thus, in order to improve scaling; fortunately, the 
wavelet-based method is insensitive to the presence of 
such zero intervals, because it only deals with the max- 
ima of T^). Due to the fact that the analyzed data is 
characterized by both the broad probability density func- 
tion and by the correlations in the temporal domain, we 
additionally estimate the contribution of each multifrac- 
tality source by calculating f(a) for the original and for 
the randomly reshuffled signals. Obviously, in the latter 
case the spectra may depend only on p.d.f.'s. 

We start with the results obtained for the DAX stocks. 
Figure 23(a) shows the fluctuation function F q (s) and the 
rescaled partition function \s'Z(q, s')| 1//<? for a randomly 
selected stock. Good-quality scaling for a broad range 
of s can be easily seen for different q's. Figure 23(b) 
and 23(c) show r(q) averaged over all 30 stocks, while 
Figure 24 displays the related mean singularity spectra 
(open circles). As this Figure documents, maxima of 
/(a) are localized approximately at a = 0.52 (MFDFA) 
and a = 0.53 (WTMM) which suggests a weak trace 
of linear correlations in data. The main difference be- 



tween the spectra obtained by each of the two methods 
is that their widths are different: MFDFA produces a 
thinner spectrum than does WTMM. We can say that 
for the stock market data WTMM detects a richer mul- 
tifractality than does MFDFA. This also refers to the 
reshuffled signals, which of course do not show any lin- 
ear correlations and their maxima are located precisely 
at a = 0.5. For the DJI stocks (Figure 26, open cir- 
cles), f(a) has a maximum at a = 0.51 (MFDFA) and 
a = 0.53 (WTMM), while for the reshuffled signals a 
maximum is at a — 0.5, which is consistent with the 
corresponding results for the German market. This con- 
sistency is also seen in the relation between widths of the 
real-data spectra derived by MFDFA and by WTMM. 
In Figure 26 both methods give spectra that on aver- 
age are narrower than for the German stocks in Figure 
24, while the opposite refers to f(a) for the reshuffled 
signals (filled squares). For DAX, it is even disputable 
whether the reshuffled-data spectra may be considered 
multifractal. However, significant standard deviations of 
data points (denoted by error bars in Figures 24 and 26) 
might question validity of these observations. Thus, we 
are not justified to infer any decisive conclusions about 
similarity or dissimilarity of the outcomes of MFDFA and 
WTMM for our financial signals. (This can be compared 
with the outcomes of an analysis of river run-off and pre- 
cipitation data in ref. |10| , where both methods produced 
results resembling each other.) 

Table 1 summarizes the results for the original and for 
the randomized signals coming from both markets. In 
each case the width Aa is larger for an original than 
for a randomized signal which is due to the existence of 
strong nonlinear correlations. A more careful inspection 
of Table 1 allows us to state that MFDFA identifies the 
temporal correlations as the principal source of multiff ac- 
tality (Aa rand < Aa). The outcomes of WTMM also go 
in this direction but they are burdened with larger sta- 
tistical uncertainity. 



V. CONCLUSIONS 

Our aim was to verify which of the two existing mod- 
ern methods of detecting multiffactality, MFDFA or 
WTMM, gives more reliable results when applied to a 
few specific sorts of data. In order to quantify their be- 
haviour we tested these methods on computer-generated 
time series with the exactly known fractal properties and 
then we performed calculations with an exemplary real 
data from a stock market. Our results indicate that from 
a global point of view the Multifractal DFA works bet- 
ter in the majority of situations presented here. First of 
all, as our examples of the fractional Brownian motion 
and the truncated-Levy process show, MFDFA is more 
reliable in properly detecting monofractal and bifractal 
behaviour than is the wavelet-based method. In such a 
case WTMM can spuriously suggest multifractality if a 
too wide range of the Renyi parameter is used. Thus, in 
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an analysis of a signal with unknown fractal properties, 
if one has to determine whether it is monofractal or mul- 
tifractal, one faces a fundamental problem of choosing 
a reliable range of q. What makes finding a solution to 
this problem difhcut is that for such a process one can- 
not assume which values of q give a correct result and 
which do not. From this perspective, MFDFA can be a 
much better option due to the fact that one may apply it 
more automatically, without paying too much attention 
to a choice of q. For obtaining reliable results it is suffi- 
cient if q are not chosen extremely large. Performance of 
both methods worsens with decreasing time series length, 
but typically the DFA-based method works better than 
WTMM also for shorter signals (of length N ~ 10 4 ). 

For actual multifractal data both MFDFA and WTMM 
are able to roughly assess f(a) at least for q > 0, while 
taking q <C into consideration is hazardous, especially 
for WTMM as our discussion on binomial cascades shows. 
A serious advantage of MFDFA over WTMM is a bet- 
ter stability of the former in respect to different P^'s; 
changing the wavelet can substantially affect the 

outcomes of WTMM (Figures 15, 17, and 19). Further- 
more, our findings indicate that WTMM performs poorly 
with signals comprising singularities of strength a > 1 
(e.g. the deterministic binomial cascades, the log-Poisson 
and the log-gamma cascades) : In this case it fails to de- 
tect the smoothest ones. This might be related to the 
fact that WTMM has inherent problems with such sig- 
nals [Mill: By construction, in order for WTMM to be 
capable of detecting the singularities of strength m m i n , 
the corresponding wavelet has to be at least ipi™™™). In 
consequence, if one is interested in investigating a broad 
range of singularities including the ones for a 3> 1, one 
has to choose a high order Gaussian derivative j/n" 1 ) or 
has to switch to MFDFA otherwise. Finally, accuracy of 
MFDFA is sensitive to the presence of linear correlations 
in data. Our results for the fractional Brownian mo- 
tion show that the more persistent the signal, the worse 
MFDFA performance. But even for H 3> 0.5 it overper- 



forms its wavelet competitor, as it does also for data with 
broad p.d.f. (Figures 12 and 15). 

However, despite these observations that are in favour 
of MFDFA, we do not disregard the wavelet-based 
method completely. For example, it can still be suc- 
cessfully used for multifractal signals provided that one 
is principally interested in fluctuations associated with 
moderate positive g's. As our analysis reveals, for such 
q's in a number of model situations WTMM acts as well 
as does MFDFA. If one is careful enough to confine an 
analysis to some small range of \q\ for which scaling of 
Z(q, s) is particularly good and based on a good statistics 
of data, then the increasing arms of the f(a) spectra for 
monofractal signals can be defined on a narrow range of 
a. Thus, if compared with some benchmark like the frac- 
tional Brownian motion, the results can indicate the char- 
acter of the signals correctly. WTMM can also be a good 
alternative to MFDFA for signals which are defined on 
a non-compact support and comprise a significant num- 
ber of zeros. Furthermore, one has to keep in mind that 
across our analysis we applied only a single wavelet fam- 
ily of the Gaussian derivatives. Although these wavelets 
are particularly popular in studies of empirical signals, 
we cannot exclude a possibility that some other type of 
wavelets can act better if applied to WTMM. Finally, a 
caution is needed about susceptibility of WTMM several 
factors in its numerical implementation. These in partic- 
ular involve noise in maxima-detecting wavelet transform 
computation at small scales and the related boundaries. 





DJI: MFDFA 


DJI: WTMM 


DAX: MFDFA 


DAX: WTMM 


Aq 


0.13 ±0.05 


0.25 ±0.10 


0.16 ±0.05 


0.26 ± 0.07 


Aa rand 


0.04 ±0.02 


0.08 ± 0.02 


0.02 ±0.01 


0.05 ± 0.02 


Aa rand /Aa 


0.30 ±0.19 


0.32 ±0.15 


0.16 ±0.07 


0.19 ±0.09 



TABLE I: Numerical results for the mean singularity spectra 
for the DAX and DJI stocks. 
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Figure 1: Example of a monofractal function: the devil's staircase obtained 
by integrating the Cantor measure, (a) Analyzed signal (top) and its wavelet 
transform (bottom) are presented together with (b) the singularity spectrum 
f(a) (main) and the scaling exponent r(q) (inset). The wavelet used in this 
calculation was ip^ . 
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Figure 2: Classical Brownian motion H = 0.5: (a) Fluctuation function F q (s) 
(open symbols) and rescaled partition function (sZ(q, s/41)) 1 / 9 (filled symbols) 
for different values of Renyi parameter: q = — 5 (squares), q = —2 (triangles up), 
q = 2 (circles), and q = 5 (triangles down). Note that the plots for different q's 
were vertically shifted in order to improve readability. Calculations were carried 
out on time series of length N = 130, 000 by applying and . Functional 
dependence of mean scaling exponent r(q) was derived from F q (s) in MFDFA 
(b) and from Z(q, s') in WTMM (c) and compared with its theoretical form 
(solid line). 
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Figure 3: Classical Brownian motion H = 0.5: Mean singularity spectra f(a) 
(open circles) obtained by using MFDFA (P( 2 > , left column) and WTMM (3) , 
right column) procedures for different range of Renyi parameter q: —3<q<3 
(top), —5 < q < 5 (middle), and —10 < q < 10 (bottom). Time series of length 
TV = 130, 000 were used; error bars indicate standard deviation of data points 
calculated from 10 independent realizations of the process. 
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Figure 4: Classical Brownian motion H = 0.5: Fraction R^ of the total number 
of signal segments (MFDFA) and of the total number of the maxima lines / e 
L(s') (WTMM) effectively used in calculation of F q {n) for MFDFA (p( 2 \ left 
column) and of Z(q,s') for WTMM (V' (3) , right column). Plots for different 
negative values of q are shown: q = —3 (top), q = —5 (middle), and q = —10 
(bottom), with three different independent process realizations of length N — 
130,000 (denoted by distinct symbols and lines). Note that for both methods 
only such ranges of n and s' that were associated with the central stable parts 
of were considered in calculation of r(q). 
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Figure 5: Classical Brownian motion H = 0.5: f(a) (open circles) for different 
time series lengths: N = 15, 000 (top), N = 65, 000 (middle), and N = 130, 000 
(bottom). Performance of MFDFA {P {2 \ left column) and WTMM (^ (3) , right 
column) is compared. Error bars indicate standard deviation of data points 
calculated from 10 independent realizations of the process. 
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Figure 6: Classical Brownian motion H = 0.5: /(a) (open circles) obtained 
with MFDFA for different polynomials: from (top left) to (bottom 
right). Time series of length N = 130, 000 were used; error bars indicate stan- 
dard deviation of data points calculated from 10 independent realizations of the 
process. 
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Figure 7: Classical Brownian motion H = 0.5: /(a) (open circles) obtained with 
WTMM for different wavelets: from ipW (top left) to tp^ (bottom right). Time 
series of length N = 130, 000 were used; error bars indicate standard deviation 
of data points calculated from 10 independent realizations of the process. 



18 



1 

0.5 


1 

0.5 


1 

0.5 


1 

0.5 




i 1 1 1 1 i 



i 1 1 1 1 



I I I I I I I I I I I I I I I I I I 



,(2) 



H — I — I — I — — I — I — I — I — — I — I — I — I — — I — I — I — h- 



J3) 



I I I I I I I I I I I I I I I I I I 



MFDFA 



i .... i 



l 1 1 1 1 l 



l 1 1 1 1 
.(1) ' 



¥ 



i i i l i i i i l i i i i l i i i i 



¥ 



(2) 



H — I — I — I — — I — I — I — I — — I — I — I — I — — I — I I I 



¥ 



(3) 



I I I I I I I I I I I I I I I I I I 



(4) 



¥ 



WTMM 



i .... i 



0.25 0.5 0.75 1 0.25 0.5 0.75 1 

a a 



Figure 8: Antipersistent fractional Brownian motion H = 0.3: /(a) (open 
circles) for different polynomials P"> of MFDFA (left column) and for different 
Gaussian derivatives of WTMM (right column). Time series of length 

TV = 130, 000 were used; error bars indicate standard deviation of data points 
calculated from 10 independent realizations of the process. 
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Figure 9: Persistent fractional Brownian motion H — 0.75: /(a) (open circles) 
for different polynomials P^ of MFDFA (left column) and for different Gaussian 
derivatives V> (m) of WTMM (right column). Time series of length N = 130,000 
were used; error bars indicate standard deviation of data points calculated from 
10 independent realizations of the process. 



20 



(a) 



1 

0.5 

£ 
^ 1 

0.5 



(b) 





1 

0.5 



5. o 

^ 1 

0.5 







H 1 1 1 1 h 



N= 15,000 

i 



MFDFA 



~i 1 1 1 r 

I 



H 1 1 1 1 H 



N= 65,000 
H = 0.3 I #=0.75 



o 



H 1 1- 



H h 



N= 15,000 

i 



WTMM 



H 1- 



^=65,000 



H = 0.3 



# = 0.75 





0.25 0.5 0.75 0.25 0.5 0.75 1 

a a 



Figure 10: Fractional Brownian motion: (a) MFDFA with and (b) WTMM 
with tp^: Mean singularity spectra (open circles) for different time series 
lengths: N = 15, 000 (top) and N = 65,000 (bottom) for the fractional Brown- 
ian motion in its antipersistent (H — 0.3, left) and persistent (H — 0.75, right) 
variants. Error bars indicate standard deviation of data points calculated from 
10 independent realizations of each process. Plots should be compared with te 
respective panels of Figures 8 and 9, where signals of N — 130, 000 were used. 
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Figure 11: Levy process qj, = 1.5: (a) Fluctuation function F q (s) (open sym- 
bols) and rescaled partition function (sZ(q, s/41)) 1 / 9 (filled symbols) for differ- 
ent values of Renyi parameter q = —5 (squares), q — — 2 (triangles up), q = 2 
(circles), and q = 5 (triangles down). Calculations were carried out on time 
series of length N = 250, 000 by applying and ijj^ ■ Functional dependence 
of mean scaling exponent r(q) was derived from F q (s) in MFDFA (b) and from 
Z(q,s') in WTMM (c) and compared with its theoretical form (solid line). 
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Figure 12: Levy process ai = 1.5: /(a) (open circles) obtained with MFDFA 
procedure (left column) for different polynomial orders: from (top) to 
(bottom) and with WTMM procedure (right column) for different wavelets: 
from i/jW (top) to -i/^ 4 ' (bottom) compared with the corresponding exact theo- 
retical values of /(a) (filled squares). Time series of length N = 250, 000 were 
used; error bars indicate standard deviation of data points calculated from 10 
independent realizations of the process. 
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Figure 13: Levy process oil — 1.5: /(a) (open circles) for different time series 
lengths: N = 100, 000 (top) and N = 250, 000 (bottom) together with the cor- 
responding theoretical values of f(a) (filled squares). Performance of MFDFA 
(P^ 2 ), left column) and WTMM (^ 3 \ right column) is compared. Error bars 
indicate standard deviation of data points calculated from 10 independent real- 
izations of the process. 
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Figure 14: Deterministic binomial cascade a = 0.65: (a) Examples of fluctua- 
tion function F q (s) and rescaled partition function (sZ(q, s/41)) 1 / 9 for different 
values of Renyi parameter: q = —5 (squares), q = —2 (triangles up), q = 2 
(circles), and q = 5 (triangles down), (b) and (c) Scaling exponents r(q) for 
MFDFA (PW) and for WTMM (ip^) are presented together with their theo- 
retical behaviour (solid line). Cascade was generated in k = 17 stages and its 
length was equal to N = 131, 072 points. 
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Figure 15: Deterministic binomial cascade: Singularity spectra (open circles) for 
different mass allocations between subintervals: a — 0.55 (top), a — 0.65 (mid- 
dle), and a = 0.75 (bottom). Left column: results from MFDFA with different 
polynomials: P^ (squares), P^ (triangles up), P^ (diamonds), and P^ 
(triangles down). Right column: results from WTMM with different wavelets: 
ipW (circles), ifj^ (triangles up), (diamonds), and (triangles down). 
Theoretical /(a) spectra are denoted by solid lines. Cascade was generated in 
k = 17 stages and its length was equal to N — 131, 072 points. 
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Figure 16: Deterministic binomial cascade with a — 0.55 (top), a = 0.65 (mid- 
dle), and a — 0.75 (bottom): Singularity spectra (open circles) for different 
cascade stages: k = 14 (N = 16,384, circles), k — 16 (N — 65,536, squares), 
and k — YJ (N = 131,072, triangles) together with the corresponding theoreti- 
cal spectrum (solid line). Left column: results from MFDFA with P^ 2 K Right 
column: results from WTMM; for each value of parameter a a wavelet with the 
best performance (see Figure 15) is used: if)^' (top), ijj^ (middle), and ip^ 
(bottom). 
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Figure 17: Log-Poisson binomial cascade with 7 = 1.4: Mean singularity spectra 
f(a) (open circles) obtained with MFDFA (left column) for different polyno- 
mials: from (top) to (bottom) and with WTMM (right column) for 
different wavelets: from ip^ (top) to (bottom), compared with the corre- 
sponding theoretical spectrum (solid lines). Cascade was generated in k = 17 
stages (N = 131,072). Error bars indicate standard deviation of data points 
calculated from 10 independent realizations of the process. 
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Figure 18: Log-Poisson binomial cascade with 7 = 1.4: /(a) (open cir- 
cles) for different cascade stages: k — 14 (N — 16,384, top), and k = 16 
(N = 65,536, bottom) together with theoretical spectrum (solid lines). Per- 
formance of MFDFA {P {2 \ left column) and WTMM (V> (3) , right column) is 
compared. Error bars indicate standard deviation of data points calculated 
from 10 independent realizations of the process. Plots should be compared with 
the corresponding panels of Figure 17, where the results for signals of k = 17 
(N = 131,072) are shown. 
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Figure 19: Log-gamma binomial cascade with 7 = 1, /3 = In 2: f(a) (open 
circles) obtained with MFDFA (left column) for different polynomials: from 
(top) to (bottom) and with WTMM (right column) for different wavelets: 
from f/jW (top) to ip^ (bottom), compared with theoretical spectra (solid lines). 
Cascade was generated in k = 17 stages (N = 131,072). Error bars indicate 
standard deviation of data points calculated from 10 independent realizations 
of the process. 
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Figure 20: Log-gamma binomial cascade with 7 = = In 2: /(a) (open 
circles) for different cascade stages: k — 14 (N = 16, 384, top) and k = 16 
(N = 65,536, bottom) together with theoretical f(a) (solid lines). Perfor- 
mance of MFDFA (P< 2 ), left column) and WTMM (^ (3) , right column) is com- 
pared. Error bars indicate standard deviation of data points calculated from 
10 independent realizations of the process. Plots should be compared with 
the corresponding panels of Figure 19, where the results for signals of k = 17 
(N = 131,072) are shown. 
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Figure 21: Log- normal binomial cascade with A = 1.1: f(a) (open circles) 
obtained with MFDFA (left column) for different polynomials: from PW (top) 
to p( 4 ) (bottom) and with WTMM (right column) for different wavelets: from 
V'*- 1 ^ (top) to ip^ (bottom) compared with the theoretical spectra (solid lines). 
Cascade was generated in k — 17 stages (N = 131,072). Error bars indicate 
standard deviation of data points calculated from 10 independent realizations 
of the process. 
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Figure 22: Log- normal binomial cascade with A = 1.1: f(a) (open circles) for 
different cascade stages: k = 14 (N — 16,384, top), and k = 16 (N = 65,536, 
bottom) together with theoretical f(a) (solid lines). Performance of MFDFA 
(P^ 2 ), left column) and WTMM (?M 3 \ right column) is compared. Error bars 
indicate standard deviation of data points calculated from 10 independent real- 
izations of the process. Plots should be compared with the corresponding panels 
of Figure 19, where the results for signals of k = 17 (N = 131, 072) are shown. 



33 




Figure 23: (a) Transaction-to-transaction logarithmic price increments of ex- 
emplary stock traded on Deutsche Borse: Volkswagen (VOW): Examples of 
fluctuation function F q (s) and rescaled partition function (sZ(q, s/41)) 1 / 9 for 
different values of Renyi parameter: q = —5 (squares), q = —2 (triangles up), 
q = 2 (circles), and q = 5 (triangles down), (b) and (c) Mean scaling exponent 
7(g) for MFDFA (P^) and for WTMM (t^ 3 )) calculated from stocks corre- 
sponding to 30 companies of DAX (empty circles). Straight lines in bottom 
panels denote a monofractal, linearly uncorrclatcd reference signal. 
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Figure 24: Mean /(a) spectra obtained by MFDFA (a) and by WTMM (b) for 
30 stocks corresponding to DAX companies (empty circles) together with their 
randomly reshuffled counterparts (full squares). Error bars for the original data 
denote standard deviation of data points calculated from 30 stocks. 
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Figure 25: (a) Transaction-to-transaction logarithmic price increments of ex- 
emplary DJI stock: IBM Corp.: Examples of fluctuation function F q (s) and 
rescaled partition function (sZ(q, s/41)) 1 / 9 for different values of Renyi param- 
eter: q = — 5 (squares), q — — 2 (triangles up), q = 2 (circles), and q = 5 
(triangles down), (b) and (c) Mean scaling exponent r(q) for MFDFA (P^ 2 )) 
and for WTMM (ip^) calculated from stocks corresponding to 30 DJI com- 
panies (empty circles). Straight lines in bottom panels denote a monofractal, 
linearly uncorrelated reference signal. 
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Figure 26: /(a) obtained by MFDFA (a) and by WTMM (b) for 30 stocks 
corresponding to DJI companies (empty circles) together with their randomly 
reshuffled counterparts (full squares). Error bars for the original data denote 
standard deviation of data points calculated from 30 stocks. 
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